Anisotropic Spectra of Acoustic Turbulence 
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We found universal anizopropic spectra of acoustic turbulence with the linear dispersion law 
ui(k) — ck within the framework of generalized kinetic equation which takes into account the finite 
time of three-wave interactions. This anisotropic spectra can assume both scale-invariant and non 
scale-invariant form. The implications for the evolution of the acoustic turbulence with non-isotropic 
pumping are discussed. The main result of the article is that the spectra of acoustic turbulence tend 
to become more isotropic. 



1. INTRODUCTION AND GENERAL 
DISCUSSION 

Wave turbulence, which describes the behavior of a 
spatially homogeneous field of random dispersive waves, 
has led to spectacular success in our understanding of 
spectral energy transfer processes in plasmas, oceans and 
planetary atmospheres [GJ. In the case of small level of 
nonlinearity e (for example, for the surface waves this is 
the ratio of the wave amplitude h to the wavelength A, 
e = h/X), there is a consistent description of the weak 
wave turbulence in terms of so called kinetic equation 
(KE) which describes the energy transfer due to inter- 
actions of three (in some cases four) waves with the con- 
servation of energy and momenta 

u(k) ± cj(ki) - u(k 2 ) = 0, fc±ki-fc 2 =0. (1) 

These equations on a classical language are a condition of 
time-space resonance. Equations (|l|) show that acoustic 
waves are special: interacting waves have to have parallel 
wave vectors: 



k II fci || k 2 



(2) 



and therefore the interacting acoustic waves are foliated 
into noninteracting rays of waves propagating in different 
directions. Immediately a few important questions arise: 

1. What are the mechanisms responsible for an energy 
redistribution between different rays? 

2. How does the energy become shared between neigh- 
boring rays? 

3. Does the energy tend to diffuse away from the ray 
with maximum energy or can it focus onto that ray? 
In the latter case, one might argue that shock for- 
mation may again become the relevant process es- 
pecially if the energy should condense on rays with 
very different directions. And finally: 

4. Is the approximation of the KE adequate for a de- 
scription of acoustic turbulence even in the case of 
small nonlinearity? 



The positive answer on the last question is problematic 
jlj . Approximation of KE is based on the randomization 
of phases of noninteracting waves leading to the Gaussian 
statistics and requires weakness of the interaction e 2 «l 
to ensure the closeness the wave statistics to Gaussian- 
ity. This happens in some physical situations, but not 
for acoustic waves. The physical reason is the above de- 
scribed foliation: for each particular (noninteracting) ray 
of waves one may pass into co-moving (with the sound 
velocity c) reference system in which all the waves are 
in the rest and therefore their interaction time goes to 
infinity. Therefore for any small level of nonlinearity the 
waves have enough time for the finally large deviation of 
phases from the random Gaussian distribution. This is 
exactly the reason why in one dimensional case one has a 
creation of shock waves (which may be described by the 
Burgers equation) for any (small) level of nonlinearity. If 
really acoustic waves tends to focus, the approximation 
of KE is problematic even qualitatively. In this paper we 
will show that this is not the case and therefore the ap- 
proximation of KE may serve at least as a starting point 
to describe non-dispersive acoustic turbulence. 

The answers on the first three questions for weakly dis- 
persive acoustic waves was done a long ago by V. L'vov 
and G. Falkovich g|. They showed that weakly dispersive 
waves, say with the dispersion law 
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really tends to focus. Namely the isotropic solution of 
the KE 



no 



(fc) = afc" 9/2 , ( a is a dimensional constant) , (4) 



found by V. Zakharov and R. Sagdeev in 1970 U is un- 
stable in the sense that the anisotropic solutions of the 
KE found in @: 
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increases with k from the pumping scale 1/L toward to 
the depth of the inertial interval kL 1 faster then the 
isotropic one (Q). In (Jfy n is the unite vector which is 
determined by anisotropy of the pumping. Observe that 
the ratio 
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n (k) 
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increases faster for more anisotropic shapes (larger I and 
p) and for more dispersive waves with bigger S. In this 
sense the non-dispersive waves with (5 = 0, w(k) oc k are 
marginal, ni p (k)/no(k) oc k° =const. and any angle dis- 
tribution n(k) is a solution of the KE at S — 0. Clearly 
this is a consequence of the foliation described above. 
However physical intuition tells as that there should be 
some mechanism of redistributing energy between neigh- 
boring rays even for fully non-dispersive waves. Indeed, 
the fact that all interacting waves have parallel wave vec- 
tors follows from the time-space resonance conditions (Q) . 
These conditions are valid with some accuracy which is 
determined by the life time of the waves and therefore 
neighboring rays can really interact. In our recent paper 



Q we generalized the KE for acoustic waves to account 
the finite width of the resonances. Generalized Kinetic 
Equation (GKE) for the "occupation numbers of waves" 
n(k, t) has the form: 



dn(k, t) 

at 



St fe ({n(fc',*)}) 



(7) 



where the collision term Stfc({n(fc', t)}) is a functional of 
the occupation numbers n{k' , t)} with all wave vectors k' 
but at the same moment of time t [which for the short- 
ness we will skip from the arguments: n(k',t) =>■ n(fe )]. 
The collision term for GKE is very similar to that for 
the KE: it is proportional to the square of the amplitude 
of three wave interactions V(k, ki, k 2 ), bilinear in n(k') 
and actually contains one three-dimensional integration: 



St,({n(fc'))) 



1 \V(k, k u k 2 )\ 2 [n(fe 1 )n(fe 2 ) - n(fc)[n(fci) + n(k 2 )]} 
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\V(k 2 , fci, k)\ 2 [n{k 2 )[n{k x ) + n(Jfe)] - n(fe)n(Jfei)] ' 



[oj(k) +uj(k 1 ) - w (fc 2 )] 5 
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(8) 



This collision term differ from that of the KE in the finite 
width rfck 1 k 2 of the resonances (pi): 



kkik2 



7(fe) + 7 (fc 1 )+7(fc 2 ) 



(9) 



where "f(k) is the damping of monochromatic wave with 
given k. This allows interaction of waves from different, 
but neighboring rays. We will see that for small nonlin- 
earity, e 2 <C 1 the characteristic angle of interaction is 
small: A6k <C 7r. This helps us to significanlty simplify 
the collision integral. 

The paper is written as follows. In Sect. || we sim- 
plify the collision term of the GKE by using the differ- 
ential in ang le approximation. Observe from Eqs. ( [22] , 
^9| ) and ( |30| ) below that the collision term St/- in the 
differential approximation has just one-dimensional Re- 
integration along the ray with direction of k and op- 
erators of differentiation in two orthogonal directions. 
Note, that Eq. ( |30| ) written for distributions n(k) with 
characteristic angle width much larger then the inter- 
action angle. This equation has an isotropic solution 
no(fc) oc fc~ 9 / 2 which coincide with the solution (|j) of 
the KE. 

In Section [| we linearize differential form of the GKE 
assuming that the deviation ni(fc) of the distribution 
n(k) from isotropic solution no(fc) is small and expand 

iii(k) = n(k, cos#) , cosd =k-n/k 

into series of Legendre polynomials fe(cos6>): 



(fc,cos0) =^/*(fc)P*(cos< 



(10) 



After that our problem foliates into set of decoupled 
equations for the Legendre polynomial with a given or- 
der I. In such a way we reduce the dimensionality of 
the problem to dimension one. Now the unknown func- 
tion fe(k) depends only on one variable k and the cor- 
responding collision terms involve only one-dimensional 
fc-integration. Our observation is that the equations for 
different Pi involve only one combination of t he param- 
eters, namely e 2 £(£ + 1). In Subsection 3C we found 
scale-invariant solutions of these equations: 



fi{k) oc 



xi> = 6 + 



ln[e 2 l(l+ l)B] 
ln(fc»L) 



(11) 



which are valid for the region of parameters where 5 < 
xi < 6. In Eq. ( |ll|) B is some number of order one, L 
is outer scale of turbulence and fc* is the wave vector for 
which 7 (A;*) is about the smallest frequency of waves in 
the system, ~ c/L. Note that the isotropic solution (^) 
no(fc) oc A: -9 / 2 and therefore when kL increases the ra- 
tio fe(k)/no(k) decrease at least as \j\fkL. It means 
that in a cascade of energy transfer from anisotropic re- 
gion of pumping down to the depth of the inertial in- 
terval energy tends to diffuse between all the rays and 
asymptotically, in the limit k — > oo acoustic turbulence 
become fully isotropical. This phenomenon of isotropisa- 
tion of non-dispersive acoustic turbulence contrasts with 
self-focusing of weakly dispersive acoustic turbulence dis- 
covered in and see Eqs. (||, ||, ^|). Note also that for 
non-dispersive waves the rate of isotropization, i.e. the 
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rate of decreasing the ration fe(k)/no(k), depends only 
on the combination e 2 £(£ + 1) and increases both with e 2 
and £. 

In Section ^ we found non-scale invariant solution of 
the linearized equations for fi{kh) which depends on 
some characteristic length A: 



fe(kA) 



\^t 2 £{£+l)C{kk) 2 



(12) 



Here C is some constant of order unity. Solution ( [l2] ) is 
valid for fcA > 1, the value of A has to be found from 
a matching of this solution with a solution for smaller 
k. Again, the rate of isotropisation depends only on the 
combination e 2 £(£ + 1) and increases both with e 2 and £. 

Note that the choice between two found solution is del- 
icate issue and depends on the various parameters of the 
problem: e 2 and £{£ + 1) separately, on the value of un- 
derground dispersion of the waves, etc. We do not think 
that it is reasonable to study this question in general, 
without referring to a specific physical realization of the 
acoustic turbulence. 

The main qualitative message of this paper is that in 
spite of anisotropic pumping at large scales L acoustic 
turbulence became to be more and more isotropic with 
increasing of its wave vector. The rate of isotropisation 
increases with increasing the level of nonlinearity e 2 and 
depends on the characteristic angle of the distribution 
A0 ~ n/£ like 1/(A6») 2 . 



2. DIFFERENTIAL APPROXIMATION OF 
THE GENERALIZED KINETIC EQUATION 

Our starting point is the generalized kinetic equation, 
GKE (0, H) which describes the interaction of neighboring 
rays due to the finite width of the three-wave resonance 
(|9|), determined by the damping increment of an indi- 
vidual wave with given k, j{k). The value of "/(k) was 
calculated in 01: 



j(k) = vk 2 , v 



A 2 N 

47TC 



(13) 



where v is the effective viscosity, N is the total number 
of the waves in the system 



N= J n{q)q 2 dq. i 
l/L 



(14) 



Integrating (g) with the help of the delta-functions, 
one has 



St fe ({n(fc')» = 



A 2 k 



(2tt) 3 



i r A^_r_ 

- / dfcifci(fc-fci) fi2 +r2 



n r 

dkiki(k + fci)^2~TW 



Here we have used the short-hand notations: 



(16) 



r± = i# 2 + fc 2 + |fc±fc 1 | 2 ], (17) 

N± = n(k ± fci)[n(fci) ± n(k)] - 7i(fc)raf>i), 

n± = c (k±k 1 - |fe±fci|) . 

Under the above assumption of dominance of the 
collincar interaction one may take ki || A; in T± to get 



r± = 2v(k 2 + kf± kfo) 



(18) 



However we cannot use the same approximation in the 
Eqs. for f2± and N±. To understand this let us have 
a look at the angular integrals in (|l^) at given value of 
k\ . It is clear that the main contribution to the integrals 
comes from the region where |f2± | varies from to a value 
of order j(k). Within this region T± is constant (18) with 



accuracy of order j(k)/ck ~ e . Similar considerations 
shows that we have to account a variation of N± within 
the region of angular integration in (fl6| ). 

To evaluate the integrals in Eq. (16) one has to estab- 



lish the relations between different wave vectors appear- 
ing in the equation. Let n be some physically unique 
direction and assume for simplicity an axial symmetry 
around n. The vectors k and k\ are almost collinear 
with a small angle 8 between them and arbitrary ori- 
ented. Then we introduce a Cartesian coordinate system 
with the axis z along k ( Fig. m). 




and A ~ yfc/p characterizes the three-wave interaction 
amplitude 



V(k, k u k 2 ) = A^Jkk x k 2 . 



(15) 



Here we accounted that the interaction is dominant by 
the interaction of almost collinear wave vectors and ne- 
glected the angular dependence of V(k, ki,k 2 ). 



FIG. 1. The coordinate system used in calculation of the 
integrals in Eq. (7). Vector n denotes the physically unique di- 
rection. The axis z is along k, the axis x is in the plane 
(k,n). The angles between vector n and vectors k and fci 
are denoted by 9 and Ox, respectively. 5 is the small angle 
between k and fei- The azimuthal angle of fci is denoted by 
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Given the geometry, we now derive a differential ap- 
proximation for the anisotropic GKE. To this end, con- 
sider (|l7|) and expand Q± and N± in S . Clearly, 



a. 



ckki 



2(k±k 1 ) 



whereas N± may be written as: 



N± = N. 



(o) 



(i) 



N. 



(2) 



Here 



N. 



(0) 



nk±q{n q ± Hfe) - UfcTlq 



(19) 



(20) 



(21) 



St M ({rc(fc')}) = ^ f ' 



dqq 2 (k-q) 2 N { ° ) (28) 
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+ / <% 2 (fc + g) 2 ivf 



Term Stfc^ ({"■(&')}) is responsible for the angular evo- 
lution (and disappears for the isotropic distributions of 
n(k), as expected): 



St fe)2 ({n(fc')}) 
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2 r 



dqq 2 {k-q) 
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2 Jo 

dqq 2 (k + q) 



N" 



(29) 



with q = fefei /fc. Terms with iV^ cx 8 cos 9 or <5 sin 6* and 

disappear after integration over <j>, while A^ 2) = 5 2 N£. 
After (free) integration over <j>, N± reads: 



3. WEAKLY ANISOTROPIC SPECTRA 
OF ACOUSTIC TURBULENCE 



(n q ± r7, fc )V 2 L n fc±fcl 



(22) 



-(nfe±q - nfc)V 2 L n fcl ± 2(Vj_nfe 1 ) ■ (V±n fc ± fcl ) 



where 







suit 



9 cos 6*i 



(23) 



and 0i is the angle between vectors fci and n; derivatives 
are taken at 6-y = 9. 

Note, that all the ^-dependence of the integrand in Eq. 
( |ilf ) is hidden in f2± and N± and therefore the integrals 
may be factorised. Using smallness of 6 we write explic- 
itly dki = Trq 2 dqdS 2 and consider the integrals over 6 2 : 



I± 



dS 1 



N±T± 



(o) 



where 



dS 2 



r(2) 



7riv£r ± y 



<5 4 (14) 2 



<5W 



r 2 ± 



r 2 ± 



r(2) 
<± ) 



2 ' 



7T N£r±c± 

2 (O'i) 2 



(24) 



(25) 



(26) 



Here the upper limit in the second integral b is deter- 
mined by the next terms of expansion of frequency in 
6, generally speaking b = 0(1). In Eq. (|26| ) C± ~ 
ln(0^/r ± ). Finally, substituting (§3] - ||) into Eq. @ 
we get the anisotropic GKE (0) with the collision term 
in the differential approximation 

St fe ({n(fe')}) - St M ({n(fe')}) + St fe>2 ({n(fe')}) . (27) 

Here Stk t o({n(k')}) originates from N±^ and coincides 
with the collision term of KE [ED : 



In this Section we will find and analyze a steady state 
weakly anisotropic solution to the anisotropic GKE (Q) 
with the collision term ( |28| , p9| ) in the differential approx- 
imation: 

St fc ,o({n(fc')}) + St fc , 2 ({n(fe')}) = . (30) 



A. Linearization of the basic equation 

Due to the weak anisotropy of the problem, we are 
seeking a solution in the form: 



n(k) — n(k, cos 9) = no(fc) + ni(fe, cos 9) 



(31) 



Here no(fe) is given by (||) and is a steady state solution 
of the isotropic problem: 



St fc ,o({rao(fc')}) = 



(32) 



The anisotropic correction is assumed to be small: 
ni(k,cos9) <C n (k). Then we substitute ( |3l"| ) to ( |30| ) 
and linearize (30) by keeping only terms proportional to 
ni and discarding terms with higher orders of the correc- 
tion. To this end, we expand the anisotropic correction 
ni(k,cos9) in a series ( |10| ) in which P^(cos#) is the Leg- 
endre polynomial of the order m, satisfying the equation: 



[Vl +£(£ + l)]P £ [cos(0)] =0. 



(33) 



Consider first St fcj0 ({n(fc')}) given by Eq. fl28|). Ac- 
cording to (^) rio(fc) is the steady state solution of the 
GKE. Therefore the terms proportional to no (ft) have to 
disappear. The result (linear in n{) is given by 

oo 

St M ({n(fc')} = J2 p dcos9)$ {k,h) , (34) 



where 
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^- ± J dqq 2 (k - q) 2 {f e (k ~ q)[n (q) - n (k)} + n (k - q)[f e (q) - f t (k)] - f e (k)n (q) (35) 
- 

CO 

- n (k)fe(q)} + J dqq 2 (k + q) 2 {f e (k + q)[n (q) + n (k)} + n (k + q)[f e (q) + f ( (k)} - ft(k)n Q (q) - n (k)f e (q)} 


To linearize Stk^({ n (k')}) we substitute in Eq. ( p9[ ) the distribution n(k) defined by (Q, [l0|,[3l]) and using ( p3[ ) to get: 

oo 

St fc , 2 ({n(fe')}) = -^^+l)P4cos(0)]$ 2 (fc,/ £ ), (36) 

l=i 

where 



,r_£ 



and 



^± £) - J{[no(g)±no(fc)]/i(fc±g)+[no(*;±g)-no(fc)]/i(?)} 



(37) 



The steady state weakly anisotropic GKE therefore 
reads: 



L 



x — 5 



00 

5^[*o(A,/i)-^ + l)*a(*,//) 



J 4 2 n o (fc)0 £ (2a; - 9) 
27rc|(a;-4)(x-5)| 



»5— a: 



(38) 



x > 5 , 
4 < x < 5 . 
a; < 4 . 



(41) 



The particular solution of this equation which may sat- 
isfy any boundary conditions may be found only if each 
term in the sum vanishes, i.e. 

$ {kJ t )=l{l+l)$ 2 {k,h) . (39) 

This is the basic equation for our study. 

B. Evaluation of the collision integral 



In order to solve Eq. (|39J) in the leading order on the 
class of scale invariant functions 



For x > 3 two integrals in ( |35j ) diverge in the IR regime 
(for q — > and for q — > k). However the leading divergent 
terms are canceled and the region of IR divergence be- 
comes x > 4. Moreover, first subleading terms also can- 
celed and real region of IR divergence in ( |35| ) is x > 5. In 
this regime one has to cutoff the integral at outer scale L, 
see first line in Eq. (fll|). For 4 < x < 5 the sum of inte- 
grals in (|35|) converge both in IR and ultraviolet (i.e. for 
5 > fc) regimes. The corresponding evaluation is given 
by the second line in Eq. (f4l|). Finally, for x < 4 the 
value of $0 is dominated by the ultraviolet (UV) diver- 
gent contribution to the second integral in ([35|) and have 
to be regularised by some UV cutoff fc* . Corresponding 



Mo) = — 



is a prefactor , 



(40) 



one has to evaluate the integrals in $0 and $2 and to 
find a leading contribution in the various regions of the 
for exponents x. These is done in Appendix |A|, here we 
only will present and discuss the results. 



evaluation is given by the last line in (41). The origin of 
fc* is different for different physical system and will not be 
discussed here. Observe from equation (^l|) that (i) $0 
has definite sign, (ii) $0 = at x=9/2, consistent with 
0, (iii) prefactor in $0 diverges for x — > 5 and x — > 4. 

In Appendix A 2 we evaluated the leading contribu- 
tions to integrals in Eq. (^) for $2 for different values 
of x. The answers may be summarized as follows: 



As we show in Appendix A 1 , the leading contributions 
to $ (pq) may be written as 



*a(*,/<) 



A 2 vCn (k)4> e 



2tt 



(42) 
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(for 6 < x) , 



for 3- < x < 6 . 



k 3 L x ~ 3 
k 3 L x - 3 + k 6 - x 
-{kLf' 2 k 6 - x + k 6 ~ x for x < ^ . 

In the IR regime we have two different contribution to 
$2, the first is divergent for any x, the second is diver- 
gent for x > 3. These contributions coincide at x = 11/2. 
One may recognize the contributions from the IR regimes 
by the presence of L in the corresponding terms. In the 
integral ( |37| ) for <E> 2 we have UV-divergence for x < 6, in 
this region the upper cutoff again is fc* and corresponding 
terms involve this parameter. Note that in contrast to 
evaluations (|df) we do not have here a window of locality 
where the collision integrals converge. 



C. Solution of the Homogeneous Equation 



Fortunately, the exponents yo = y% = § coincide in 
the window 5 < x < 6, meaning that at least the k- 
dependence of &o(k,fe) and <j> 2 (fc, /^) is the same. Next 
promising observation is the coincidence of sings: both 
function are positive. It means that we have a chance 
to satisfy Eq. ( |39| ) by a proper choice of the exponent x. 
Indeed, according to ([id], |I|), Eq. ( |39| ) for these values of 
x becomes: 



L x ~ 



B' (x-5)l(£ + l) 



vkl 



(44) 



where B' is some positive, i- independent, dimensionlcss 
factor [we believe of O(l)] which accumulated all un- 
known factors in our estimates. 

Let us estimate effective viscosity v. Substituting 
isotropic solution (|J) in ( fL4"| ) one gets N ~ aL 3 / 2 . Now 
from (|l^) one has 



aL 3 ' 2 /p 



(45) 



Our goal now is to find the solution of the homoge- 
neous Eq. (|39|). We seek the solution in the scale invari- 
ant sector. Denote the scaling exponents of &o(k, fi) and 
§2{k 1 fe) as yo and y 2 : 



$o(fc, ft) cc 4>tk- ya , $ 2 (fc, f t ) oc 4>ik- 



112 



(43) 



The x dependences of these exponents are determined by 
Eqs. d| H) and shown in Fig. |. 

$0 converge for 4 < x < 5 and diverge for other values 
of x (in IR for 5 < x and in UV for x < 4). $2 never 
converge with different rate of divergence in IR and UV 
limits (two branches, see Fig. ||). The ratio between IR 
and UV-divergent terms in $ 2 depends on the wave vec- 
tor k, therefore, generally speaking we have to account 
both IR and UV terms. 
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Next one has to relate the dimensional factor a with the 
dimensionless amplitude of waves e. To do this we eval- 
uate the total energy of the acoustic waves 



d 3 ku{k)n (k) 



Air ca 

l/L 



dk 
W 2 



caVI . (46) 



We define e 2 as the the ratio of the total energy in the 
acoustic-wave system E &c to the total kinetic energy of 
the media -Ekin which is about p c 2 . Parameter e may be 
treated as the ratio of the amplitude of velocity in acous- 
tic waves to the mean square velocity in the media caused 
by its kinetic energy in the thermodynamic equilibrium. 

Now we have e 2 ~ a-Jhjpc which together with ( |45| ) 
gives 



e 2 cL 



Using this evaluation in Eq. ( |44| ) one has: 
L x - 6 = B" <l 2 i{£ + l)kt~ x , 



(47) 



(48) 



where B" absorb one more unknown factor from the 
evaluation of v. Therefore the solution is achieved for 
x = xq,i, where 



xo.e = 6 + 



ln[e 2 ^ + 1)B"] 
ln(fc* L) 



(49) 



FIG. 2. Scaling exponents yo (solid line) and j/2 (dashed 
lines) for different values of x. Arrows show the directions of 
divergence for the corresponding integrals. "IR" and "UV" 
stay for the infrared and ultraviolet limits, "locality" denotes 
the region of convergence of $0. 



These exponents give the fc-dependence of the solution 
( p9| ) in the depth of the inertial range. However, the 
functions (j)£ remain unknown. To find them one have to 
match the inertial range solution to the pumping at the 
IR boundary (k = l/L). 



G 



D. Effect of Inhomogeneous Terms 

Sometimes the inhomogeneous terms in the KE causes 
additional solutions which may play important role in 
the evolution of spectra in k. Consider first the origin of 
the inhomogeneous terms. Linearizing Eq. ( |34| ) we have 
concluded that n§ contribution to St^.o — is zero. This 
is true only if the IR limit of the integral is indeed zero, 
which is not the case. In the energy containing interval 
< k < L^ 1 we have a non-universal behavior of n(fe) 
which has to be accounted. A simple way to do this is 
first to evaluate the contribution of this region as 



Stfe i0 ,inhom({"-(fc'}) 



A 2 no(fc) 

^-'i 

2irc\/L 



(50) 



Note that in this region n(fc) is not isotropic and thus 
Stfc ] o,inhom({"-(fc'}) must depend on the direction of k. 
Expanding this dependence into the spherical harmonics 
(in our case of the axial symmetry of the problem into 
the Legendre polynomials) we have an inhomogeneous 
contribution to the Eq. (p9): 



in horn 



(k,n ,. 



A 2 n (k)n , e 
2ttcVL 



oc k z 



(51) 



where non-universal numbers no/ may be related to 
via "boundary conditions" at fc « X/L; 



n 



TX-9/2i 



(52) 



An important observation is that the scaling exponent 
z of &i n hom(k, no/) in (51) is independent of x (because 
this term is independent of fi(k) altogether). Moreover, 
z = i, which is exactly the same as exponents yo and 
y2 of the homogeneous part of the GKE 4>o and <&2 ■ It 
means that for any k one has to account the contribution 



of $inhom(fc, fto/) m the balance equation (p9). Repeat- 
ing the same calculations as in deriving Eq. (j44[) with the 
help of Ens. (|3S| , fll| ) and (^) and accounting now eval- 
uations ( pl| , |52|) we are approaching again Eq. (44) with 
additional contribution to the factor B"; their sum may 
be denoted as B. Physically it means that resonant (i.e. 
with the same scaling exponent) inhomogeneous term has 
shifted the exponents xg relative to their "homogeneous 
values" xq^ (|49|). By replacing the unknown constants 
B" => B we obtained from (^) final equation ([ll|) for 
the scaling exponents of xg . 



4. NON-SCALE INVARIANT SOLUTIONS 

In previous Sect, we found all scale- invariant solutions 
and observe that they have finite region of applicability 
for which their scaling exponents 5 < xg < 6. Outside 
of this region we have to look for non-scale invariant so- 
lutions which contain explicitly a characteristic length- 
scale A which will be chosen later. To do that we have 
to choose fi(kh) in some reasonable form, for example 



ft(n) = a ^K x 



1 



+ 



"4,£ 



(53) 



exp [ — bin yl \ 



kA> 1. 



where x, y, a,/ and bg are some unknown numbers. Ob- 
viously there is no UV divergences in $o and $2 for 
such choice of fg. Furthermore, since in leading order 
lim fg = do there is no IR divergence in $0 and $2 asso- 

ciated with fg, so that the IR behavior of the <&o and $2 
is dominated by no(g , IR ), q m ~ 1/L. 

Then in expression ( |35| ) we keep only the terms pro- 
portional to no(g) and ng(fc — q): 



$ (fc,/^)~A_[i j dqq 2 (k-q) 2 {fg(k-q)n a (q)+n (k-q)[fg(q)-fg(k)}-fg(k)n (q)} (54) 
+ J dqq 2 (k + q) 2 {fg(k + q)n Q (q)-fg(k)n (q)}' . 
Now let us change variables in the first integral q — ► k — q. Then we can evaluate integrals at q = q 1B only: 



A 2 

.2 



4irc 



i J dqq 2 (k - q) 2 {fe{k - q)n (q) + n (k - q)[fg(q) - fg(k)} - fg{k)n (q) + fg{q)n (k - q) 

+n Q (q)[fg(k -q)- fg(k)} - fg(k - q)n Q (k -<?)}+ f dqq 2 {k + q) 2 [fg{k + q)n Q (q) - fg(k)n (q)} 

A 2 k 2 8 2 fg(k) 



4ttc 



dqq 2 k 2 2n (q)[fg(k + q) - 2f e (k) + fg(k - q)} 



8ttc dk 2 



dqq n (q) . 



(55) 
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The second derivative of fi is the result of the usual dou- 
ble cancellation of the IR divergence in the isotropical 
part of KE. 

Now we have to evaluate $2 at q <C k for the choice of 
ft given by (p3[). Let us first evaluate N" for q <C k: 



N'!_ {1) ~ \{n Q {q)h{k-q) + [n Q {k-q)-n Q {k)]h{q)}, 

(56) 



5. SUMMARY 



Let us summarize the logic of this paper: 

• We start from Generalized Kinetic Equation 
(GKE) §), 



N'L {1) ~ \n (q)h(k), g«fc 

Similarly, 

N'f> ~ \n Q {q)f t {k), 
Now use (|Aq ) to estimate $2: 
A 2 



®2(k,fe 



16tt 2 c 
A 2 



2 2 Avk 2 £ 

dqq k n Q {q)fi{k) 

cq 



16tt 2 c IR 



q 2 kHe 2 Ln (q IR )f e (k), 



(57) 
(58) 

(59) 
1/L, 



where we used estimation (Efl) for effective viscosity v. 
Actually we have not calculated accurately the q — > fc 
contribution in $2- Note however, that this contribu- 
tion would be of lower (in q) order, because lim^fc ~ 
-[c/c 2 /2(fc-g]l. 

Now Eq. ( p9[ ) acquires the form (introducing dimen- 
sionless constant C, generally speaking of the order of 
unity) 



Q^k 2 n (q ln ) 



d 2 h(k) 
dk 2 



Ce 2 L£(£+l)q 2 k A n (q 1B )f t (k) 



(60) 



Note that no(q m ) cancels from both sides of this equa- 
tion. Further notice that both terms are of positive sign. 
Canceling q 1B and k from both sides and substituting 
k = fc/A with q 1B ~ 1/L we get 



d 2 fi(n) 

dK 2 



C£(£+ l)e 2 K 2 fi(n) , K=kL 



(61) 



where we chose A ~ L. This equation may be solved in 
the Bessel functions (of order j of imaginary argument). 
For our goals, however, it would be enough to present 
only its asymptotic form for k> 1: 



0,2,'. 



+ 



exp 



(62) 



where coefficients 02/, 04, t, etc. may be found iteratively: 



«2 



e = 3/y/Ce(£ + l)( 



but these terms really are be- 



yond the accuracy of Eq. ( |6l| ) itself. Solution (|62j) shows 
that anisotropic corrections decay exponentially, the rate 
of decay increases with e 2 and £(£ + 1). 



• The conditions of time-space resonances (Q) dic- 
tates almost collinear propagation, so we expand 

in the transverse direction, and arrive to the 
differential approximation for the collision term in 
GKE @ with @ and @; 

• We assume weak anisotropy with anisotropic cor- 
rection in the factorized form (|3l|); 



• After substitution ( pl| ) in Eqs. ( p7| , [2q , p9| ), we 2m- 
earize resulting three-dimensional equations (which 
are differential in angles and integral along the 
rays); 



• Expansion of these equations in series in Legendre 
polynomials Pi foliates the 3-dimensional problem 
in decoupled sets on 1-dimcnsional ones, for each £ 
separately; 



• We have found the power-law solutions of this 
equations (jll]) and the exponential solutions ( p^ ) 
which are governed by the same parameter e 2 £{£ + 
1). Thscs solutions describe the phenomenon of 
isotropization of non-dispersive acoustic turbulence 
in the cascade process of energy transfer from 
anisotropic region of pumping down to the smaller 
and smaller scales (fc — > 00). In the limit k — > 
00) the statistics of acoustic turbulence approaches 
isotropical flux equilibrium. 



The main result of the article is that the spectra of acous- 
tic turbulence tend to become more isotropic. 
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APPENDIX A: ANALYSIS OF CONVERGES 
OF INTEGRALS IN THE COLLISION TERMS 

1. <&o(fc, fe) behavior 

a. IR-regime 

In Eq. ([35]) one has three different regions dangerous 
for IR divergence. These are the regions q — > and q — » k 
in the first integral and q — > in the second one. The first 
integral ^ Jq ■ ■ ■ may be split into two integrals i J^ 2 . . . 
and ^ J fe , 2 ... In the second one we may change dummy 
variable q —> q' = (k ~ q) and then to re-denote q' — > (7. 
After that one sees that the integral | J,, 2 . . . equal to 

fc/2 , ,1 r ,1 n , • , 1 • /T7P1\ 1 r*; 



l/o • • • and therefore the first integral in (|35|), | J Q . . .. 

may be rewritten as J* .... This integral in the region 
q <C k together with the second integral in (35) may be 
written as 

*(>(*:, ft) « ^ y fc/ ^ 2 (fc - <z) 2 /£(?) (ai) 

x [n (fc - g) - n (g)] + (fc + g) 2 /*(s)M* + g) - n(fc)] . 

Substituting fe(q) from (|I]) one observes here a usual 
double cancellation of the leading and the first sublead- 
ing terms in the IR regime. Therefore integral ( |Al[ ) may 
be evaluated (up to a factor) as follows: 



A 2 n Q (k) 



A-/2 



dqq 4 fi{q) 



(A2) 



Now one sees that the integral ( |A2| ) converges for x < 5 
and diverges for x > 5. In the later case we take the 
outer scale L as an IR cutoff. Then it is easy to see that 
$0 may be evaluated in different regions of x as: 



A 2 n Q {k)(j) l 

2lTC 



k 5 - x /(5-x) (for x < 5), 

ln(fci) (at a; = 5) , (A3) 

L x ~ 5 /(x - 5) (for 5 < x) . 



Note, that in this limit and for x close to 5, $0 > 0. 



b. UV-regime 

Next we analyze the region q^> k, i.e. the UV-regime 
in the second integral in Eq. (|35|). Here the most dan- 
gerous terms reduce to 



*o(*,/<) 



^4 2 no(fc) 

4-7TC 



dqq 4 [fi(k + q) - f e (q)] . (A4) 



There is a cancellation of the leading contribution in the 
region q 3> k. Therefore this integral converges at x > 4. 
In a compact form the UV-evaluation of <E>o may be writ- 
ten as: 



$0 



A 2 kn (k)4>£ 
2nc\x - 41 



x > 4 . 
x < 4 



(A5) 



Here fc* is an UV cutoff of the integral which may have 
different nature for different physical situations and will 
not be further clarified here. Note, that in the UV-regime 
for x close to 4, $0 < 0. It is known also, that $0 = 
for x — 9/2, the scaling index o f the isotropic solution of 
the KE. Now we can join Eqs. (A3) and ( |A5| ) and write 
evaluation of the leading contribution to $ = in all 
regions of x in the form (41). 



2. $2(fe, fe) behavior 

a. IR-regime 

Unlike &o(k, fe) , there is no double cancellation of the 
IR contribution (i.e. in the limits q — > and q — > k) in 
the sum of the integrals ( p7| ) for $2(^1,/^). The reason 
is that Q'± is not invariant under q — > k — q transforma- 
tion. Now, in the IR region q -C k one has the following 
simplifications: 



r± 


= 2v(k 2 + 


k\ ± fejfei) ~ 


n± 


^5 2 n'i, 




c± 







2z//c 2 , 
ckk\ 

z ± fe] 

= £ 



(A6) 



±cg/2 



Substituting these equations into (|3j]) one gets in the IR 
region: 



®2(k,fe 



vk 4 A 2 C 
2tt 2 c 2 



_ k 



dqq 



- N" W - N 

2 " 



(A7) 



The most dangerous terms in the combinations iV^ are: 

A^f } ~ i{no(g)//(A ± q) + [n (k ± q) - n (k)]f e (q)} ■ 

(A8) 



Together with Eq. ( |32[ ) these yields 
j{no(q)fi(k) T ||n 



at 



//(f) 



(*)/<(«)} 



(A9) 



Then in the IR limit $2 reads 
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vk"A 2 L 



_ fe 



(A10) 



A 2 W2 

27T 2 C 2 fc 



dqq 6 N" W 



-n (q)fi(k) + -j-n (k)ft(q) 

where we did not care about the numerical factor, just 
carrying the signs. Since no(k) oc fc~ 9 / 2 the no(q)fe{k) 
term always diverges; the corresponding contribution to 
the integral behaves as (kL) 5 ^ 2 k 2 no(k)fg(k). The second 
term will diverge if fi(q) — 4>e/q x with x < 3. Symboli- 
cally, its contribution to the integral are: 



no(k) 4>e 



L' 



(for x < 3) , 
(for x > 3) . 



(All) 



Out of N'!^ the most dangerous for UV term is 



N'l {l) ~ - [n (fc)(/,(fc + «) - ft(q))] ~ --r-n (k)f e (q) , 



4q 



(A15) 



where we used that fe(q) oc q~ x and q ^> k. Now in the 
UV regime the $ 2 (fc, /^) term may be evaluated (up to a 
numerical factor) as 



Observe that for x > 11/2 the contribution of the second 
term in (AlC) dominates the contribution of the first one. 
Finally, the IR contributions to $2(^7 fe) are summarized 
as: 



Mk,M 



vk 3 A 2 £n (k)(f) e 



2tt 2 c 2 

-(kL) r °/ 2 k 3 - x ( for x < 11/2) 



(A12) 



-x-3 



(for x > 11/2) 



6. UV-regime 

Next we have to consider the UV regime of the integral 
in ( |37| ) i.e. the region of integration with q —> kf k. 
Here we have the simplifications: 



r+ = 2v{k 2 + k\ ± fcfci) ~ 2i/g 2 , 

cfcfci 



(A13) 



Now $2 becomes 



n'i = ± 



2(k±ki) 
= t . 



ck 

T 



A 2 
8tt 2 c 



\^vq i 



dqq 2 (k + q) 2 ^r 



-N' 



(A14) 



®2{k,ft 



A 2 vCn (k) 



dqq 5 fi{q) 



(A16) 



This integral converges for x > 6 and may be written as 
follows: 



A 2 j/£n o (/c)0£ 



fcj x for a; < 6 , 



2^ 2 c 2 |x-6| fc6 _x for x>6 



(A17) 



Combining Eqs. ( A17 ) and ( A12 ), one has after some 
minor manipulations (neglecting factors of order one, dif- 
ference between C and C, etc.) Eq. (142). 
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